Relativistic BCS-BEC Crossover at Zero Temperature 
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We investigate the BCS-BEC crossover at zero temperature in the frame of a relativistic model. 
The universality of the BCS-BEC crossover for non-relativistic systems breaks down in relativistic 
case and the crossover can be induced by changing the density. When the effective scattering length 
is much less than the fermion Compton wavelength, we recover the non-relativistic result if the gas 
is initially in non-relativistic state. At ultra-strong coupling where the scattering length is of the 
order of the Compton wavelength, a new BEC state appears. In this state the condensed bosons 
become nearly massless and anti-fermions are excited. The behavior of the Goldstone mode and 
the mixing between the amplitude and phase modes are significantly different in different condensed 
regions. 

PACS numbers: ll.30.Qc, 12.38.Lg, 11.10. Wx, 25.75.Nq 



I. INTRODUCTION 

It was found many years ago that the Bardeen-Cooper- 
Shriffer (BCS) superfluidity/superconductivity in degen- 
erate Fermi gas and the Bose-Einstein condensation 
(BEC) of composite molecules can be smoothly con- 
nected when Eagles and Leggett found that the ground 
state wave functions of BCS and BEC states are essen- 
tially the sameP, Q. The BCS-BEC crossover phenom- 
ena can be realized in high temperature superconduc- 
tivity and atomic Fermi gas where the s-wave scattering 
length can be adjusted [3,lIL IE @|- 

Recently, the study on BCS-BEC crossover has been 
extended to relativistic Fermi systems with the Nozieres- 
Schmitt-Rink (NSR) theory @, and the Bose- Fermi 
model[9[. In the NSR theory above the critical tem- 
perature, a new feature of BCS-BEC crossover is 
found: There exist two kinds of BEC states[7|, the 
non-relativistic BEC (NBEC) and the relativistic BEC 
(RBEC). In the RBEC state, the condensed bosons be- 
come nearly massless, anti-matter can be excited, and 
matter and anti-matter are nearly degenerate. A nat- 
ural question is: Does the RBEC appear only in rela- 
tivistic systems and never exist in non-relativistic sys- 
tems? In this paper, we will argue that the RBEC phase 
can happen in any system if the attractive coupling is 
large enough, even though the system is initially in non- 
relativistic state. From text books ^fj|. when the Fermi 
momentum of a system is much smaller than the fermion 
mass, we can treat the system non-relativistically. This 
statement is true only if the coupling is weak enough. 
For sufficiently strong coupling, the system will become 
relativistic even though the Fermi momentum is much 
smaller than the fermion mass. 

We will extend in this paper the theory of BCS-BEC 
crossover in symmetry breaking state [5[ at zero temper- 
ature to relativistic systems. The requirement for such 
an extension is that we should recover the non-relativistic 
theory in a proper limit. We will focus on two relativistic 
effects: the anti-fermion degrees of freedom and the non- 
trivial fermion mass. The first effect leads to the appear- 



ance of the RBEC state at sufficiently strong coupling, 
and the second effect breaks the universality of BCS-BEC 
crossover. We will also investigate the collective modes 
evolution in the BCS-BEC crossover. 



II. RELATIVISTIC THEORY OF BCS-BEC 
CROSSOVER 

The physical motivation why we need a relativistic the- 
ory of BCS-BEC crossover is mostly due to the study of 
QCD phase diagram, especially the dense quark matter 
which may exist in compact stars and can be created in 
heavy ion collisions. However, we will point out that the 
theory is also necessary for non-relativistic systems when 
the coupling is strong enough. To this end, let us first 
review the non-relativistic theory of BCS-BEC crossover 
in a dilute Fermi gas. 

Leggett 's mean field theory [l| is successful to describe 
the non-relativistic BCS-BEC crossover at zero temper- 
ature. For a dilute Fermi gas with fixed density n = 
kj/(3ir 2 ), where kf is the Fermi momentum, the BCS- 
BEC crossover can be found if one self-consistently solves 
the gap and number equations for the pairing gap Ao and 
the fermion chemical potential 
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where E k and £k are defined as E^ = 
£k = k 2 /(2m)-/i n , and a s is the s-wave scattering length. 
The fermion mass m plays a trivial role here, since the 
BCS-BEC crossover depends only on the dimensionless 
parameter r\ — 1/ (kfa s ), which is the so-called universal- 
ity for non-relativistic syetems. The BCS-BEC crossover 
can be characterized by the chemical potential \x n : It 
coincides with the Fermi energy e/ = fe?/(2m) in the 
BCS limit rj — > — oo but becomes negative in the BEC 
region. In the BEC limit rj — > +oo, one has ji n — > —Eb/2 
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with E\y = \/(ma 2 s ) = 2rfef being the molecular bind- 
ing energy. Therefore, in the non-relativistic theory the 
chemical potential will tend to be negatively infinity in 
the BEC limit. 

A problem arises if we look into the physics of the 
BEC limit from a relativistic point of view. In the 
relativistic description, the fermion dispersion becomes 



! ± n r , where =F correspond to fermion and 



anti-fermion degrees of freedom, and fi r is the chemical 
potential in relativistic theory ll| . In the non-relativistic 
limit |k| -C to, if fi r is of the order of m, we can neglect 
the anti-fermion degrees of freedom and recover the non- 
relativistic dispersion £^ ~ k 2 /(2m) — (fi r — m). There- 
fore, not \i r itself but the quantity [i r — m plays the role 
of non-relativistic chemical potential While the 

chemical potential fi n can be arbitrarily negative in non- 
relativistic theory, \x r is under some physical constraint 
in relativistic theory. Since the molecule binding energy 
can not be larger than two times the constituent mass, 
the absolute value of the non-relativistic chemical poten- 
tial [i n = /x r — to = —E(,/ 2 at strong enough coupling 
can not exceed the fermion mass to, and the relativistic 
chemical potential [i r should be always positive. If the 
number density n satisfies kf <C to, the non-relativistic 
theory works well when the coupling is not very strong 
and the binding energy satisfies E\, <C 2m. However, if 
the coupling is strong enough to ensure Ej, ~ 2to, rela- 
tivistic effects will appear. From E^ ~ 1/(too 2 ), we can 
roughly estimate that the non-relativistic theory becomes 
unphysical for the s-wave scattering length a s ~ 1/m. 
This can be understood if we consider 1/m as the Comp- 
ton wavelength A c of a particle. 

What will then happen at rj — > +oo in an attractive 
Fermi gas? At sufficiently strong coupling with E\, — > 2to 



and fi r 



0, the dispersions £ k and £7 for fcrmions 



and anti-fcrmions become nearly degenerate, and non- 
relativistic limit can not be reached even for systems with 
kf <C to. This means that the anti-fermion pairs can be 
excited and become nearly degenerate with the fermion 
pairs, and the condensed bosons and anti-bosons become 
nearly massless. Without any model dependent calcu- 
lation, we can confirm an important relativistic effect 
in BCS-BEC crossover: There should exist a relativis- 
tic BEC(RBEC) stateQ which is smoothly connected to 
the non-relativistic BEC(NBEC) state. The RBEC is not 
a specific phenomenon for relativistic fermion systems, it 
should appear in any fermion system if the attractive 
coupling becomes strong enough, even though the initial 
non-interacting gas satisfies kf <C to. 

We now start to construct a general relativistic model, 
which we expect to recover the non-relativistic theory 
in a proper limit. We consider only fermion degrees of 
freedom in the original Lagrangian. The most general 
form of the Lagrangian density can be expressed as 



C = i/> (i^dp -m)tp + Ci, 



(2) 



where ip, ip denote the Dirac fields with mass to, and Cj 
describes the attractive interaction among the fermions. 



It is shown that the dominant interaction is in the scalar 
channel J p = + [H], [l3[ . For the sake of simplicity, we 
consider the zero range interaction and take the form 



(3) 



where g is the attractive coupling constant, and C — 
*7o72 is the charge conjugation matrix. Generally, by in- 
creasing the attractive coupling, the crossover from con- 
densation of spin-zero Cooper pairs at weak coupling 
to the Bose-Einstein condensation of bound bosons at 
strong coupling can be realized. 

We start our calculation from the partition function in 
imaginary time formalism, 



Z= D^D^ dT S ^ c +^^\ 



(4) 



where is the inverse temperature, — 1/T, and fi is the 
chemical potential corresponding to the net charge den- 
sity tp^tp which is determined by the charge conservation. 
Similar to the method in the study of superconductivity, 
we introduce the Nambu-Gorkov spinors 



* = ( V Vc ) 



(•5) 



and the auxiliary pair field A(x) = gip T (x)Ci'j5i['( x )/2, 
and then perform the Hubbard-Stratonovich transforma- 
tion, the partition function can be written as 



Z = / L>* D^D AD A" 



(6) 



in terms of the inverse Nambu-Gorkov propagator 

G _1 = i-y^d^ - to + ^7o(73 + «75 Afi + + ij 5 A*cr_ , (7) 

where a± = (a± ± ia2)/2 are defined in Nambu-Gorkov 
space with <Ji(i = 1, 2, 3) being the Pauli matrices. Inte- 
grating out the fermion fields, we finally derive the par- 
tition function 



Z = DAD A* e~ s ^ f 



with the bosonic effective action 



S eff - 
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(9) 



While in general case one should include the contri- 
bution from pair fluctuations [H, 0, Q, the Leggett mean 
field theory is already good to describe the BCS-BEC 
crossover at zero temperature. This can be seen from the 
fact that, the dominant contribution of fluctuations to 
the thermodynamic potential at low temperature is from 
the massless Goldstone modes and is proportional to T 4 . 
At zero temperature it vanishes [5(- While the mean field 
theory can describe well the BCS-BEC crossover, it can 
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not predict precise values of the universal constants in 
the unitary limit and the boson-boson scattering length 
in the BEC limit. Since our goal in this paper is to in- 
vestigate the BCS-BEC crossover only, we will take the 
mean field approximation in the following. 

In the mean field approximation, we consider the uni- 
form and static saddle point A(x) = Ao which satisfies 
the stationary condition SS e ff [AqI/SAq = 0. The ther- 
modynamic potential density Q — S e ff[A ]/{f3V) at the 
saddle point can be evaluated as 
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where we have defined the notations E^ = 



^ 2 + A 2 



and £^ = ek ± M with £k 



\/k 2 + TTl; 

get the gap equation for the condensate Aq 



Minimizing fi, we 
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and the number equation for the fermion density n 
k}/(3n 2 ), 
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The first and second terms on the right hand side of 
Equations (fTTj) and (| 1 2[) correspond to fermion and anti- 
fermion degrees of freedom, respectively. 

The model is non-renormalizable and a proper regular- 
ization is needed. We subtract the vacuum contribution 
dfl/dAQ\T= n =A =o from the gap equation, namely we 
replace the bare coupling g by a renormalized coupling 
U, 
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To recover the non-relativistic theory correctly, the chem- 
ical potential in vacuum should be chosen as /i = m 
rather than fi = 0. Such a subtraction is consistent 
with the formula derived from the two body scattering 
matrix[8j]. The effective s-wave scattering length a s can 
be defined as U = 4ira s /m. While this is a natural exten- 
sion of the regularization in non-relativistic theory, the 
ultraviolet divergence can not be completely removed, 
and a momentum cutoff A still exists in the theory. The 
relativistic Fermi energy Ef in the theory can be define 

i 2 , which recovers the Fermi kinetic 



as Ef — ^jk 2 



energy Ef 



ef = k 2 1 (2m) in non-relativistic limit. 



III. RELATIVISTIC EFFECTS IN BCS-BEC 
CROSSOVER 

In the non-relativistic theory, there are only two 
characteristic lengths kj 1 and a s , and the BCS-BEC 



crossover shows the universality: After a proper scaling, 
all physical quantities depend only on the dimensionless 
coupling 1] = X/(kfa s ). Especially, in the unitary limit 
a s — > oo, all scaled physical quantities become universal 
constants. Unlike the non-relativistic theory where the 
fermion mass m plays a trivial role, in the relativistic the- 
ory a new length scale, namely the Compton wavelength 
A c = m _1 appears. As a consequence, the BCS-BEC 
crossover should depend on not only the dimensionless 
coupling r), but also the quantity £ = kf/m = fc/A c . 
Since the cutoff A is needed, there is also a A/m depen- 
dence. By scaling all energies by e/ and momenta by 
kf, the gap and number equations (fTTj) and (fT2|) become 
dimensionless, 
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with the definitions E„ 



h (AoA/) 2 , & = 

e x ±fi/ef, e x = 2C 1 ^x 2 + (~ 2 and z = C _1 A/m. It be- 
comes now clear that the BCS-BEC crossover is charac- 
terized by three dimensionless parameters, 77, Q and A/m. 

We now study when we can recover the non-relativistic 
theory in the limit ( < 1. Expanding e x in powers of 
C, e x = x 2 + 2Q~ 2 + 0(C 2 ), we can recover the non- 
relativistic version of £ x . However, we can not simply 
neglect the terms corresponding to anti-fermions, namely 
the second terms on the right hand side of Equations 
(fT4l . Such terms can be neglected only when \fi — m\ 
and Ao are both much less than m. When the cou- 
pling is so strong to ensure fj, — > 0, the contribution 
from anti-fermions becomes significant. In this case 
the condition ( < 1 is not sufficient to describe the 
non-relativistic limit, and the other important condition 
should be a s <C 1/m which guarantees the molecule bind- 
ing energy E^ <^ 2m. With the dimensionless coupling 77, 
we have rj = \/{kfa s ) = (m/kf)(l/ma s ) m/kf = C ■ 
Therefore, the complete condition for the non-relativistic 
limit can be expressed as C <C 1 and ())< 1. 

To confirm the above statement we solve the gap and 
number equations numerically. In Fig[T]we show the con- 
densate Ao and non-relativistic chemical potential \i — m 
as functions of 77 in the region — 1 < r\ < 1 for several 
values of £ < 1. In this region the cutoff dependence 
is weak and can be neglected. For sufficiently small £, 
we really recover the Leggett result. With increasing £, 
however, the universality is broken and the deviation be- 
comes more and more remarkable. On the other hand, 
when we increase the coupling rj, especially for 77 — * 
the difference between our calculation at any fixed £ and 
the Leggett result becomes significant due to relativis- 
tic effects. This means that even in the case C "C 1 we 
can not recover the non-relativistic result when the cou- 
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pling 77 is of the order of This can be seen clearly 
from the 77 dependence of A and /1 in a wider 77 region, 
shown in Fig|2] We found in our numerical calculation 
with the cutoff A/m = 10 a critical coupling i] c ~ 2£ _1 
which is consistent with the above estimation. Beyond 
this point the chemical potential fj, approaches zero and 
the condensate Ao becomes of the order of the relativis- 
tic Fermi energy Ef or the fermion mass m rapidly. In 
the region around r] c , the relativistic effect becomes sig- 
nificant, even though the initial non-interacting gas is in 
a non-relativistic state. 
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FIG. 1: The condensate Ao and non-relativistic chemical po- 
tential n — m, scaled by the non-relativistic Fermi energy e/, 
as functions of 77 in the region — 1 < r\ < 1 for several values 
of £. In the calculations we set A/m = 10. 



We can derive an analytical expression for the critical 
coupling 7] c or U c for the RBEC state. At the critical 
coupling, we can take fi ~ and Ao <C m approximately, 
and the gap equation becomes 
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from which we have 
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FIG. 2: The condensate Ao and chemical potential /x, scaled 
by the relativistic Fermi energy Ef, as functions of 77 in a wide 
region — 1 < r\ < 20 for several values of £. In the calculations 
we set A/m — 10. 



with the definitions f(x) = ln(x + \/x 2 + 1) and Uq = 
4ir/m 2 — AttX 2 .. While in the non-relativistic region with 
77 -C C the solution is almost cutoff independent, in 
the RBEC region with 77 ~ the solution becomes 
sensitive to the cutoff. 

To see what happens in the region with 77^ ~ 1, we 
discuss the fermion and anti-fermion momentum distri- 
butions 7i_(k) and n + (k), 



n±(k) 



(17) 



In the non-relativistic BCS and BEC regions with 77C <C 
1, we have Ao <C Ef, and the anti-fermion degrees of free- 
dom can be safely neglected, 7i + (k) ~ 0. In the BCS limit 
with Ao <C £/, 7i_ (k) deviates slightly from the stan- 
dard Fermi distribution at the Fermi surface, especially 
we have n_(0) ~ 1. In the deep BEC region, we have 



(15) A ~ r)ef and |/i - 



77 ef and in turn |/x — m\ 3> Ac 



and ti_(0) <C 1, and 7i_(k) becomes very smooth. How- 
ever, in the RBEC region with 77^^ 1, fi approaches zero, 
and the anti-fermions become nearly degenerate with the 
fermions, 
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where n±(0) can be large enough as long as Ao is of the 
order of to. 

In the non-relativistic BCS and BEC regions, the total 
net density n = n_ — n + is approximately n ~ n_ = 
J G? 3 k/(27r) 3 n_ (k), and the contribution from the anti- 
fermions can be neglected, n+ = J d 3 k/(27r) 3 n+(k) ~ 0. 
However, when we approach to the RBEC region, the 
contributions from fermions and anti-fermions are almost 
equally important. In the RBEC region with Ao < m we 
can estimate 

A ° f A ,u fc2 

A^A Am A\ . . 

= tV i- -r arctan - ■ 19 

8ir 2 \ A m / 

For m <C A, the second term in the bracket can be omit- 
ted, n_ ~ n + ~ AqA/(87t 2 ). The densities of fermions 
and anti-fermions are both much larger than n, and their 
difference produces a conserved net density. 
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FIG. 3: The phase diagram in the Uo/U — kf/m plane. The 
line is defined as fi = m. 



V. COLLECTIVE MODE EVOLUTION 



IV. DENSITY EFFECT IN BCS-BEC 
CROSSOVER 

In the non-relativistic BCS-BEC crossover, the cou- 
pling strength and number density are reflected in the 
theory in a compact way through the dimensionless quan- 
tity 77 = l/(kta s ), and changing the density of the system 
can not induce a BCS-BEC crossover. However, if the 
universality is broken, there would exist an extra density 
dependence which may induce a BCS-BEC crossover. In 
non-relativistic Fermi gas, the breaking of the universal- 
ity can be induced by a finite range interaction 1J| . In 
the relativistic theory, the universality is naturally bro- 
ken by the C = kf/m dependence which leads to the extra 
density effect. 

To study this phenomenon, we calculate the phase di- 
agram in the Uo/U — kf/m plane where Uo/U reflects 
the pure coupling constant effect and kf/m reflects the 
pure density effect. The reason why we do not present the 
phase diagram in the 77— £ plane is that both 7/ = l/(fc/a s ) 
and £ = kf/m include the density effect. The phase di- 
agram is shown in Figj3] The BEC state is below the 
line [i = m and the BCS-like state is above the line. We 
see clearly two ways to realize the BCS-BEC crossover, 
by changing the coupling constant at a fixed density and 
changing the density at a fixed coupling. Note that we 
only plot the line which separates the two regions with 
(i > m and (i < m. Above and close to the line there 
should exist a crossover region, like the phase diagram in 

The density induced BCS-BEC crossover can be re- 
alized in dense QCD such as QCD at finite isospin 
density (l5l. [Til IT7I Il8j and two color QCD at finite baryon 
density. The new feature in QCD is that the quark mass 
m decreases with increasing density due to chiral symme- 
try restoration at finite density, which would lower the 
line in the phase diagram Fig[3J 



To investigate the Gaussian fluctuations above the sad- 
dle point Ao, we write A(x) = Ao + <fi(x) and expand the 
action to the second order in 6 to obtain 
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M n (Q) = - + -y2Tr[i 75 g l 1 (K + Q)i 75 G^(K)] 



where $ is defined as $t(<?) = (<£*(Q),<£(-Q)), Q = 
(q, iv m ) with v m = 2m7r//3 (to = 0, ±1,±2, •••) is 
the four momentum of the collective mode, and ~YliQ = 
T Ylm I ^ 3( l/ (2 71 ") 3 denotes integration over the three mo- 
mentum q and summation over the frequency v m . The 
inverse propagator M of the collective mode is a 2 x 2 
matrix with the elements 

1 1 x ^ 

+ - 

9 

M 12 (Q) = i^Tr[7 75 S 12 (^ + Q)*75So 12 W] , 

1 K 

M 21 (Q) = M 12 (-Q), 

M 22 (Q) = Mn(-Q), (21) 

where K = (k, iu> n ) with u> n = (2n + l)7r//3 (n = 
0, 1, 2, ■ • • ) is the fermion four momentum, and Qq {i, j = 
1,2) are the elements of the fermion propagator Qo = 
G[Aq] in the Nambu-Gorkov space given by 
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with the energy projectors 
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At zero temperature, Mu and M12 can be evaluated as 
1 
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with 7± = 1/2 ± (k • q + e k )/(2ek£k+q) and the coherent 



coefficients «) 2 = (l+£/E£)/2 and ( vf) 2 = (1 - 
g/E£)/2. 

Taking the analytical continuation iv m — > w -M0 + , the 
dispersion u) = w(q) of the collective mode is obtained 
by solving the equation 



det M[q, w(q)] = 0. 



(25) 



To make the result more physical, we decompose the com- 
plex fluctuation field <j>(x) into its amplitude mode X(x) 
and phase mode 9(x), 4>(x) = (A(x) + i9(x))/^/2. Then 
the fluctuation part of the effective action takes the form 



A* 9* 



M+ + M12 iM u 
-iM n M^i - Mi 2 



A 



(26) 



where (q, to) = (Mu (q, ui) ± Mu (q, —uo))/2 are even 
and odd functions of lo. From M{^(q, 0) = 0, the ampli- 
tude and phase modes decouple exacrtly at lo = 0. Fur- 
thermore, M+(0,0) = Mi 2 (0,0) ensures that the phase 
mode at q = is gapless, i.e., the Goldstone mode. 

We now determine the velocity of the Goldstone mode, 
w(q) = c|q| for lo, |q| <C mink-i-E^ }. For this purpose, we 
make a small q and lo expansion of the effective action, 



M++M12 = A + C\q\ 2 -Du 
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The Goldstone mode velocity now reads 

r 2 Q 

B 2 /A + R' 



(28) 



and the corresponding eigenvector of M is (A, 9) = 
(— ic\q\B/A, 1 ), which is a pure phase mode at q = 
but has an admixture of the amplitude mode controlled 
by B at finite q. The explicit form of A, B, R and Q 
can be calculated as 
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(29) 



In the non-relativistic limit with ( < 1 and £77 <C f, 
we have |/x — m| <C m, all the terms that include anti- 
fcrmion energy can be neglected, the fermion dispersions 
E^ and £ k can be approximated by E^ and £k, and the 
function T+ can be approximated by 1. In this limit the 
functions Mu and M12 are the same as the ones obtained 
in non-relativistic theory[5]. In this case, we have 



B 
R 

Q 



1 
2 
1 

8 
I 

16 



d 3 k g k 
{2*YEl 
d 3 k 1 

d 3 k 1 



m 



A|k^ 

E 2 . m 2 



(30) 



(27) 



In the BCS limit, all the integrated functions peak near 
the Fermi surface, we have B — and c 2 = Q/R. Work- 
ing out the integrals we recover the well known result 
c = CA/3 for non-relativistic BCS superfluidity. In the 
BEC region, the Fermi surface does not exist and B 
becomes non-zero. An explicit calculation shows that 
c = Q/^/^Trq <C C0- This result can be rewritten as 
c 2 = AitnBCLB /m% where mu = 2m, as = 2a s and 
riB = tl/2 are corresponding mass, scattering length and 
density of a dilute Bose gas, which recovers the result of 
the Bogoliubov theory of dilute Bose gas with short-range 
repulsive interaction[19(. 

In the RBEC region with £77 ~ 1, we have fi — > 0, 
and the terms include anti-fermion energy become nearly 



7 



degenerate with the fermion terms and hence can not be 
neglected. Since the quantity B vanishes at fi — > 0, we 
have c 2 = limA^ooQ/-R hi the RBEC region. Taking 
fi = which leads to 



R 
Q 



i r d 3 k i 



4 7 (27r)3£3 

i r d 3 k i 

8 J J^fEj 



3- 



k 2 



A 2 k 2 



(31) 



where E'k = ^/e 2 . + A 2 , is now the degenerate dispersions 
in the limit /i — * 0, we find c = 1 in this region. This is 
consistent with the Goldstone boson velocity calculated 
from the relativistic boson field theory of BECplj. 

In the ultra-relativistic BCS state with kj 3> m and 
Aq <C /i — Ef, all the terms that include anti-fermion 
energy can be neglected again. In this case we have B = 
and 



R 



dk 



16t 2 Jo [(k - fi) 2 + A 2 ] 3/2 ' 



32tt 2 



dk- 



A 2 



o [(k - fi) 2 + A 2 



21 5/2' 



(32) 



A simple algebra shows that Q/R = 3, and hence we 
obtain the well-known result c = 1/V3 for the ultra- 
relativistic BCS superfluidity. 

The mixing of amplitude and phase modes is quite dif- 
ferent in different regions. In the weak coupling BCS 
region, all the integrated functions peak near the Fermi 
surface and we have B = due to the particle-hole sym- 
metry, and hence the amplitude and phase modes decou- 
ple exactly. In the NBEC phase with (r/ -^i 1, while the 
anti-fermion term in B can be neglected, we have B ^ 
since the particle-hole symmetry is lost, and the mixing 
is strong. In the RBEC region, while both particle-hole 
and anti-particle-anti-hole symmetries are lost, they can- 
cel each other and we have again B = 0. This can be seen 
from the fact that for fi — ► the first and second terms 
in B cancel to each other. Thus in the RBEC region, the 



amplitude and the phase modes decouple again. This 
can also be explained in the frame of the boson field the- 
ory of BEC. In non- relativistic theory, the off-diagonal 
elements of the inverse propagator are proportional to 
iuj 20], which induces a strong mixing between ampli- 
tude and phase modes. However, in relativistic theory, 
the off-diagonal elements are proportional to i^o;[Tll[2]|. 
which vanishes when BEC of massless bosons happens. 



VI. SUMMARY 

We have investigated the BCS-BEC crossover at zero 
temperature in a relativistic Fermi gas model. Unlike the 
non-relativistic theory, the fermion mass plays a nontriv- 
ial role and serves as a new length scale. As a conse- 
quence, the universality breaks down and there exists an 
extra number density effect on the BCS-BEC crossover. 
When the effective scattering length is much less than 
the fermion Compton wavelength and the Fermi momen- 
tum is much less than the fermion mass, we can recover 
the non-relativistic theory. At the ultra-strong coupling 
where the effective scattering length is of the order of 
the Compton wavelength, the RBEC state appears. In 
this state the condensed boson becomes nearly massless 
and anti-fermions are excited. The sound velocity of the 
Goldstone mode and the mixing between the amplitude 
and phase modes are quite different in different regions, 
and the results agree well with the boson field theory. 

In this paper, we investigated only the relativistic BCS- 
BEC crossover at zero temperature where all pairs get 
condensed. At finite temperature, zero momentum con- 
densed pairs can be thermally excited and one should go 
beyond the mean field theory to treat properly the non- 
condensed pairsQ. The study on the relativistic BCS- 
BEC crossover in the symmetry breaking phase at finite 
temperature is in progress. 
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